Dynamical chiral symmetry breaking in sliding nanotubes 
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We discovered in simulations of sliding coaxial nanotubes an unanticipated example of dynamical 
symmetry breaking taking place at the nanoscale. While both nanotubes are perfectly left-right 
symmetric and nonchiral, a nonzero angular momentum of phonon origin appears spontaneously 
at a series of critical sliding velocities, in correspondence with large peaks of the sliding friction. 
The non- linear equations governing this phenomenon resemble the rotational instability of a forced 
string. However, several new elements, exquisitely "nano" appear here, with the crucial involvement 
of Umklapp and of sliding nanofriction. 
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A popular high-school physics demo is a clamped os- 
cillating rope or guitar string. While being imparted a 
strictly planar vibration at one end, the string initially 
vibrates, as expected, within the plane. Above a certain 
amplitude however the plane of vibration spontaneously 
and unexpectedly begins to turn around, right or left 
with equal probability [1]. Since all along the string's 
Lagrangian (including the external forcing) remains com- 
pletely left-right symmetric (i.e., nonchiral) this is a pro- 
totype example of what may be called dynamical chiral- 
ity breaking taking place in the macroscopic world. Dy- 
namical spontaneous symmetry breaking of chiral sym- 
metry is actually rife in nature, including examples such 
as the Taylor- Couette instability in hydrodynamics 
the spontaneous chirality of vibrating cilia in biology [. 
and many other macroscopic scale examples. 

We discovered, in simulations of the frictional sliding of 
carbon nanotubes, a nanoscale example of dynamical chi- 
ral symmetry breaking in particular in two coaxial nan- 
otubes, one forced to slide inside the other. While both 
nanotubes are perfectly left-right symmetric and nonchi- 
ral, angular momentum strikingly jumps to nonzero val- 
ues in correspondence to some sliding velocities, coinci- 
dent with large phonon-related peaks of the sliding fric- 
tion. The theory of this phenomenon yields non-linear 
equations that differ from the string problem by newer 
elements that are exquisitely nano, now involving Umk- 
lapp processes and sliding nanofriction. 

We conducted classical molecular dynamics simula- 
tions of an inner, infinite and termination- free, (5,5) arm- 
chair single-wall carbon nanotube sliding at speed v in- 
side a coaxial (10,10) nanotube, as in Fig. [T^. Standard 
empirical potentials were used for the intratube [1] and 
intertube [5[ interactions. Details of simulation are given 
m Ref. 0. Temperature was controlled, usually at T = 
300 K, by an algorithm designed to preserve angular mo- 
mentum [7;]. As noted earlier by Tangney et al. [8], the 
tube-tube sliding friction is not hydro dynamical, in fact 
not even monotonic with v, but develops sharp peaks 
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FIG. 1: (a) Coaxial sliding of (5,5)@(10,10) nanotubes. (b) 
Sliding friction force per inner tube atom and angular mo- 
mentum Jcm from simulations versus speed v at T = 300 
K. Note the frictional peaks and the threshold near 780 m/s. 
The nonzero angular momentum signals nanoscale dynamical 
chiral symmetry breaking in correspondence with the thresh- 
old, and with peaks at 570 m/s and 450 m/s. A, B, and 
C designate related resonance regions, described in text, (c) 
Outer tube radial motion Fourier spectra with n — 2 (green) , 
3 (blue), and 5 (red) symmetry, showing resonant enhance- 
ment in correspondence with peaks and threshold. 



and onsets at selected speeds. We found sharp frictional 
peaks near v — 450 m/s, 570 m/s, 720 m/s, and an im- 
portant threshold onset near 780 m/s, shown in Fig. [TJd. 
These peaks are known to generally arise out of para- 
metric excitation of outer nanotube "breathing" phonon 
modes, classified by an angular momentum index n (for 
tangential quantization around the tube axis). Our peak 
positions differ from earlier ones [8] due to our lack of 
tube terminations (implying mode uniformity along the 
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FIG. 2: Tube cross-sections, showing a shape- rot at ion 
(pseudo-rotation) and an overall rotation (marked by red 
atoms) for v = 820 m/s (washboard period a/v = 300 fs). 
(a)-(g) Snapshots at intervals of 50 fs comprise one period of 
pseudo-rotation, where the pentagon shape rotates clockwise, 
(h) An overall counter rotation of nearly tt in 0.5 ns. 



tube), and also to our different inter-tube potentials. As 
noted earlier [6], the non-monotonic F-v characteristics 
implies a "negative differential friction" , whereby an in- 
creasing applied force F yields an innner-tube velocity 
that grows by jumps and plateaus, rather than smoothly. 

The surprise comes from analysing, at the frictional 
peaks, the two parts of angular momentum J around 
the tube axis y: its center-of-mass (rigid body rotation 
with angular velocity uj) and shape-rotation ("pseudo- 
rotational") parts J = Jqm + ^pseudo, with Jcm = 

Y^i m i[n x x ri)] y , and J pse udo = ^Zi m i(^i x ^i)y 
Generally zero at generic v due to lack of nanotube chi- 
rality, J pS eudo = -Jcm = 0, J pse udo jumps to nonzero val- 
ues at the frictional peaks and past the threshold, where 
jpseudo = -^cm 7^ (total J is clearly conserved), see 
Fig. [TJd. Simulations at lower temperatures (not shown) 
reveal other smaller frictional peaks, but the main find- 
ings of the present paper, including chirality breaking, 
are still present, and indeed even stronger. 

Cross-section snapshots at intervals of 50 fs, at v = 820 
m/s of Fig. [2] show a large pentagonal distortion of the 
outer tube, corresponding to an n = 5 breathing mode. 
A vertex of the pentagon moves clockwise (or counter- 
clockwise with 50-50 probability in different simulations) 
with angular velocity £l w /5, where Q w = 2i\vja is the 
"washboard" frequency associated with the sliding speed 
v and a = 2.46 A is the intertube potential periodicity. 
The pentagon rotation signals a non- vanishing J pse udo > 
0. Since the total J = 0, the center-of-mass acquires an 
equal and opposite counter-rotation Jcm = —jpseudo < 
0, as demonstrated by the red atom marker in Fig. O 

Theory - The mechanism underlying symmetry break- 
ing is the tube non-linearity, reminiscent of the string 
instability [1]. The non-linear motion of a string forced 
to vibrate along a transverse direction x at frequency Q 



is described by 

X = -[cul T + K(\X\ 2 + \Y\ 2 )]X + Fcos(Qt) 

Y = -[< T + iv(|X| 2 + |F| 2 )]r, (1) 



where Wk,T — k^/ fi/M is the bare transverse frequency 
with wavevector k (/i is the Lame constant and M the 
mass density) , X, Y are two orthogonal linear polariza- 
tions, and F is the forcing amplitude. The nonlinear- 
ity K > accounts for the increase of the transverse 
frequency at large amplitudes, where the larger over- 
all elongation leads to an effective string tension in- 
crease. Due to K > there is a critical frequency 

Q cr = ^JujI t + (KF 2 /16) 1 / 3 beyond which the string, 

although forced along x, spontaneously develops both 
X and Y modes, an elliptic polarization resulting from 
a purely linear forcing, with a spontaneous dynamical 
chirality breaking. The critical amplitude of the X- 
mode beyond which the elliptical polarization sets in 
is (2F/K) 1 / 3 (this will be demonstrated in Fig. [3)3 be- 
low). In the nanotube case, the sliding inner tube ex- 
cites phonon modes of the outer one, uniformly along 
the tube in our termination-free case. The outer (10,10) 
nanotube has, among others, doubly- degenerate modes 
with angular momentum n = ±2, ±3, ±4 and ±5, etc. 
Let u x (0) and u z (0) be displacements in the x (tangen- 
tial) and z (radial) directions at angular position 9 on the 
outer tube circumference, and u UyXy ± and u UyZy ± their re- 
spective n-mode amplitudes, travelling clockwise (+) and 
counter-clockwise (— ). Suzuura and Ando (SA) [9[ de- 
scribed these modes at the quadratic level in a contin- 
uum model. Third- and fourth-order energy terms in u 
mix the different n-modes in all possible ways compati- 
ble with conservation of angular momentum. The main 



nonlinear terms oc J dx(^r 



-)( 



du z 
dx 



%0 2 to third 



order, and oc J dx(^- — ^) 4 to fourth order turn out to 
be as crucial as in the string. 

"n = 5" excitation - The 10-chain atomic structure of 
the sliding (5,5) inner tube stimulates at all speeds v 
the outer tube with a deformation corresponding to a 
linearly polarized, k y = 0, n = 5 mode with a wash- 
board frequency Q w = 2itv/a. Away from resonance, 
< ^5 ? the excitation amplitude of the n = 5 outer 
tube mode is small, and linear. Near and above res- 
onance, however, the amplitude grows, and nonlinear- 
ities become dominant. Starting from the SA model 
[9[ including non-linearities, after some approximations 
(mainly u Xjrii ± ~ (i/n)u z ^ n ^± as appropriate from lin- 
ear eigenvector analysis, neglecting small mode-mixing 
third-order terms while keeping large fourth-order contri- 
butions which generate terms such as |^ 2)5)+ /_| 2 ii 2 ,5,+), 
we arrive at equations for the n = ±5 phonon amplitudes 
of the suggestive form 

X = -[^ 2 + Q + (P/4)(|X| 2 + |F| 2 )]X + 4Fcos(^t) 
Y = -K 2 -Q + (P/4)(|X| 2 + |F| 2 )]r, (2) 
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FIG. 3: (a) Simulation results (dots) in the threshold region 
v > 780 m/s, compared with the model (thick solid and dot- 
ted lines) of Eqs. J2}. X,Y = \u z ,5,+ =b ^* ?5; _ | are amplitudes 
of X(t) and Y(t) defined in the text. Y becomes nonzero to- 
gether with nonzero Jcm = — Jpseudo above the critical thresh- 
old speed. The hysteresis found in the model is absent in the 
simulation, (b) Forced vibrating string from Eqs. (pQ), with 
uj kjT = 1, K = 5, F = 0.002, and M 100. a, b are the 
X, y-mode amplitudes. The critical frequency Q cr = 1.0054 
and the threshold X-mode amplitude a cr = 0.093 (inset). 

where the two variables X(t) = (1^,5,+ + u* 5 _^ e tQwt + 
c.c. and Y(t) = i(u z ^ : + — u* z 5 _^e lViwt + c.c. parame- 
terize the radial displacement u z ^(9,t) in terms of two 
orthogonal standing waves, u z ^(6,t) = X(t)cos(56) + 
y(t)sin(50), and the fourth-order term P again shifts 
upwards the frequency as the amplitude increases. The 
important difference between Eqs. ([2]) and those of the 
string is the "nano" Q-term, representing an Umklapp 
process coupling two outer tube modes with n = ±5 
through a "reciprocal lattice vector" of the inner tube, 
with its 10 carbon double-chains. The Umklapp term 
splits the X/Y frequencies to approximately \J uj\ =L Q. 
A static double-tube phonon calculation confirms split 
frequencies uo x = 3.35 and uo y = 3.15 THz which measure 
the Umklapp strength. 

In nanotube sliding, unlike a guitar string that can 
be pinched soft or hard, the only controllable param- 
eter is the speed v. Nevertheless, due to the similar 
fourth-order non-linear effects, the physics of nanotube 
sliding, in the region with v > 780 m/s (shaded and 
labelled "5" in Fig. [TJd), resembles that of the string, 
see Fig. [3^. Near the 780 m/s threshold, Q w « 3.17 
THz, just above oj y = 3.15 THz, a spontaneous symme- 
try breaking occurs, the F-amplitude growing from zero, 
and the center-of-mass angular momentum with it. In 
the approximations considered so far, this should cor- 
respond to a pure "n = 5" excitation; and indeed this 
is close to reality, as shown in the Fourier spectrum of 
Fig. [Tfc, where at v = 820 m/s the most important peak 
appears at the washboard frequency Q w . Note that the 
effective n = 5 mode frequency is dragged along by the 
washboard, growing with v and with the friction magni- 



tude. When the speed grows larger than 830 m/s, simu- 
lations show the additional excitation of n = 2 or n = 3 
modes, due to third-order non-linearites not accounted 
for in Eqs. ([2]). Here the dynamics, still chiral, becomes 
more complex than that of the pure "n = 5" mode. In 
Fig. [3^ solutions of Eqs. ([2|) are shown as shaded thick 
lines, with parameters P = (1728/50)a/Mi? 4 , Q, and 
F which were fit to the bulk modulus of the nanotube 
(a = 15000 A 2 /ps 2 ), the intrinsic phonon frequencies 
and the splitting for n = 5 (Q = 50.5/ps 2 ), as well as 
the energy corrugation (F = 5 A 2 /ps 2 ) of the inter-tube 
potential. The model agrees fairly well with hard simula- 
tion results, including a crossing of X and Y amplitudes 
at around v = 815 m/s, which only occurs due to Umk- 
lapp. A feature predicted by the model is a hysteretic 
behaviour with respect to increasing and decreasing v. 
The hysteresis is not seen in simulations both at 300 K 
and at 50 K. We suspect that frictional Joule heating 
could be sufficient to remove it. Also, our model being 
mean- field in character, it does not allow for fluctuations. 
The simulations of course do, and indeed occasional re- 
versals of the sign of Jcm are observed when |Jcm| is 
very small; these reversals are suppressed, most likely by 
inertia, when |Jcm| is large. 

"2 + 3 = 5 " excitation - Consider now the main fric- 
tional peak labelled "2 + 3" in Fig. [TJd, around v « 570 
m/s. This corresponds, as shown by Fig. Sj to a mixed- 
resonance involving three modes, n = 2, 3, and 5 simul- 
taneously (see Fourier spectrum in Fig. [Tfc). The inset 
in Fig. 2] shows that the cross-section snapshot of both 
tubes in this region of v, where frictional excitation is 
huge, is mostly elliptical, indicating a dominant n = 2 
deformation (see also the large n = 2 peak in the Fourier 
spectrum of Fig. []J), with a slight triangular (n = 3) 
deformation. Our anharmonic continuum model nicely 
explains these results, the Umklapp terms (present in 
nanotubes but not in the string) playing a crucial role 
[101 ]. In this region of v we can truncate the full hier- 
archy of coupled equations to 6 equations involving the 
relevant Fourier modes u x / z ^ n with n = 2, 3, 5 [10[. To 
give a flavor of the full theory [lOj], we report here the 
equation of motion for mode u z ^, obtained under the 
linear approximation u x $ « (i/5)u Z £ [Tol ] . 

Uz,5 = -7^,5 - Qu* Zl 5 - (^5 + B rn\u z ,m\ 2 )u Zj5 

ra=2,3,5 

-A dir u z , 2 u z , 3 + A umkl < >2 < 3 + 2Fcos(Sl w t), (3) 

where the last term represents the washboard forcing 
(due to inner tube sliding), while the Q-term is the pre- 
viously discussed Umklapp term —5^5. (Similar equa- 
tions hold for the modes n = 2,3, except that the Q-term 
and the explicit forcing are missing.) The physics of the 
£? m -terms is an effective increase of the bare frequency 
due to the excitation of modes m = 2, 3, 5, via fourth- 
order non-linearities. A dir and A umkl represent the effect 
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FIG. 4: Simulation results (dots) and comparison with the 
model (with a velocity backshift of 78 m/s) for the 2+3 peak 
(thick solid and dotted lines). The model results are based 
on the 6 equations, generalization of Eqs. (|3]), for u x , n and 
u Zj n with n — 2, 3 and 5. X n ,Y n — \u ZlTl ,+ ± ^, n ,-| are 
amplitudes of X n (t) and Y n (t). Y n becomes nonzero together 
with nonzero J pS eudo above the critical sliding speed, for "2 + 
3" excitations. Hysteresis is not seen in the simulation. 

of third-order non-linearities, whereby a combined exci- 
tation of modes 2 and 3 can lead, by angular momentum 
conservation, to a term influencing mode 5 either in a di- 
rect way, via 2 + 3 = 5, or through Umklapp, via — 2 — 3 = 
—5 — > 5. The stationary solutions of the nonlinear 
coupled equations for 1^/2,71=2,3,5? obtained numerically, 
reveal a rather simple structure, namely, u x / z ^(t) = 

and u x / Zi s(t) = u x / z ^ + e lQ3t , with ^ 2 + ^3 = ^w'- in 
words, mode 5 oscillates with the driving washboard fre- 
quency Q w , while modes 2 and 3 "feel" an indirect driv- 
ing from mode 5 and oscillate at frequencies Q2 and Q3 
(renormalized by tube-tube coupling and anharmonicity) 
in such a way as to realize a perfect "resonance" with the 
driving frequency (see Fourier spectra Fig. [Tfc). 

As simulations show, the amplitudes of X n — u z , n ,+ + 
u* n _ and Y n = u Zirii + — u* n _ (in qualitative agreement 
between model and simulation) are equal to each other 
for both n = 2 and n = 3, implying u z ^/3,- = 0. That 
means, in full agreement with the model, a full chiral- 
ity for both modes, with waves travelling totally clock- 
wise (or counter-clockwise, with equal probability). The 
n = 5 X and Y amplitudes are, on the contrary, unequal 
and much smaller than n = 2 and 3. The continuum 
model actually predicts - in striking agreement with the 
simulation - that the n = 5 mode rotates opposite to 
the chirality of modes 2 and 3 (i.e., counter-clockwise, 
|w 25 5 5 _| > |i^ ? 5 5+ |, for the case shown here), an effect 
entirely due to Umklapp terms [To| . 

So far we have explained the two main features labelled 
A in Fig. [2d. The features labelled B appear to represent 
a replica, where single outer tube n = 3, 5 excitations re- 
place the double tube ones n = 3, 5. For example, in the 



spectrum for 455 m/s (Fig. [Tfc), the n — 3s mode is at 
1.14 THz, the outer tube value, different from the joint 
n = 3 mode where both tubes vibrate together. Anal- 
ogously, the n = 5s mode is excited at 720 m/s. Here, 
however, the friction is so small that the excited angu- 
lar momentum is essentially invisible. The two regions 
labelled C, finally, are irrelevant to the present context 
and will be discussed elsewhere [lOj . 

In summary we have described a spontaneous break- 
ing of chiral symmetry occurring in nanoscale friction. 
Our very large speeds and termination-free conditions 
differ strongly from experimental conditions where nan- 
otube sliding was observed so far. [11]. Attempts at ob- 
serving this chirality breaking could nonetheless be of 
considerable interest, considering in addition the pos- 
sibility to examine the effect of electronic frictional 
dragging [12] an avenue which seems worth considering. 
Although the modes that play a major role in electron- 
phonon coupling[13], and nanotube superconductivity 
[14j are different, breathing modes excitation was re- 
cently demonstrated by STM tips [15]. The physical un- 
derstanding obtained for our microscopic nonlinear dy- 
namical system might also serve as a prototype for phe- 
nomena in the field of nanomotors. Examples of sponta- 
neous chiral symmetry breaking that are relevant to the 
dynamics of nanomotors have been described in rotax- 
anes [16], where a series of chemical reactions can give 
rise to a biased Brownian motion of a small ring around 
a larger one in either direction. Our example represents 
a first idealized case where the origin of nanoscale chiral 
symmetry breaking is strictly physical. 
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